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Abstract — 

The Bayesian Land Surface Temperature estimator pre- 
viously developed has been extended to include the effects 
of imperfectly known gain and offset calibration errors. It 
is possible to treat both gain and offset as nuisance param- 
eters and, by integrating over an uninformitave range for 
their magnitudes, eliminate the dependence of surface tem- 
perature and emissivity estimates upon the exact calibration 
error. 
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Sea Surface Temperature. 
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I. Introduction 

AS a practical matter, land surface temperature (LST) 
estimates retrieved from radiances reported by a re- 
mote sensor will be subject to some unavoidable level of 
calibration error, which need not be very accurately known. 
While all forms of quantitative exploitation of radiometric 
data in remote sensing are afflicted to some degree by cali- 
bration error effects, this problem may be a special concern 
for the Bayesian multiband LST algorithm £Q . That is be- 
cause the Bayesian algorithm iterates on a range of plausi- 
ble surface temperatures, within which the estimated LST 
value is obtained as an expected value. Should uncom- 
pensated calibration errors lead to a temperature interval 
which does not bracket the true surface temperature, the 
algorithm in its present form has no way to recover, and 
may return a surface temperature estimate with degraded 
accuracy. 

This note sketches the extension of the Bayesian ap- 
proach to LST retrieval to include effects of a simple form of 
unknown calibration error. After a review of the Bayesian 
approach to LST retrieval, the calibration error is param- 
eterized as linear in the true aperture radiance. The joint 
prior probalility for the calibration error parameters is then 
obtained by imposing the requirement that two distinct 
observers agree on its mathematical form. Finally, the 
Bayesian LST estimator is extended to include linear cal- 
ibration error by treating the calibration error parameters 
as nuisance variables, and integrating them out of the final 
estimators for surface temperature and emissivity. 

II. Elements of Bayesian LST estimator 

The Bayesian land surface temperature retrieval algo- 
rithm is developed in 1 , which may be consulted for de- 
tails. The approach to LST retrieval presented in that 
earlier paper consists of three elements: 



f . The forward model for sensor aperture radiance, as- 
sumed linear in surface emissivity: 

I F (k) = ek B k (T)exp(~-) + F^ (0)exp(——) + /J(t, /i) 

(1) 

If. (t, /x) and (0) are the upwelling diffuse radiance at 
nadir optical depth r (top of the atmosphere, or TOA, 
for spaceborne sensors; fi is the cosine of the angle with 
respect to zenith) and the downwelling irradiance at the 
surface, respectively. B k {T) is the Planck function at sur- 
face temperature T. The emissivity is e k , and the surface 
reflectance p k — 1 — e&. Note that QJ assumes Khirchoff's 
law; this is done solely for simplicity. It is also assumed (at 
least initially) that the sensor has high spectral resolution. 

2. The MAXENT form of the conditional probability 
of observing radiance I 0,0,01 m the presence of noise- 
equivalent radiance a: 



P(I\T,e k ,a) 



exp 



2ct 2 



dl 
a 



(2) 



3. The prior probability of surface temperature and 
emissivity pQ: 



P(T,e k \K) = ^dTde k 



(3) 



The posterior probability for the surface temperature 
and emissivity, given observed radiance and available 
knowledge, is obtained from these quantities by use of 
Bayes' theorem: 

P{T, e k | I,K) oc P{T, e k \ K)P{I \ T, e k , a). (4) 

III. Calibration error model 

By hypothesis, the physical radiance Iq at wavenumbcr 
k is related linearly to the reported radiance J: 

7o = (l + a)/ + /3 (5) 

Both a and p are assumed to be small quantities; 



and 



a « 1 



P « h 



(6) 



(7) 
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It is the physical radiance Jo which goes into 0J. Recalling 
that the forward model Q is linear in e, the exponent in 
that expression is 

(I - If) 2 = (((1 + a)I + 13)- (Ae + B)f (8) 

which is quadratic in a, (3, and e. 

IV. Prior probability for calibration error 

PARAMETERS 

In order to obtain a useable estimator, it is necessary 
to find the prior probability for the calibration error pa- 
rameters a and f3 0,0,0- As in 0, two equally cogent 
observers must relate their descriptions of radiance, and of 
calibration error, by a Lorentz transformation [H|,[H1 con- 
necting one (primed) coordinate description moving with 
velocity v along the observation axis with respect to the 
other (unprimed) one by 



k = r]k 



(9) 



where the Doppler factor r\ is given in terms of the boost 
parameter 

' (10) 
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by 



- (^) 2 

i] = 7(1 — v/c) 

j l-v/c 
~ V 1 + v/c 



(11) 



The quantity rj is real and nonvanishing for physical 
Lorentz transformations. Let 



P(a,/3 I K) = g(a,(3)dad(3 



(12) 



be the prior probability assigned by Vladimir in the un- 
primed frame, and 



P(a',/3' I K) = h(a',(3')da'dp' 



(13) 



be that assigned by Estragon, viewing radiance in the 
primed frame. The prior probabilities in the two frames 
are related by 



h(a',(3')da'df3' = ,r 1 g{a 1 f3)dad(3 



where 



J = det 



d(a',(3') 



(14) 
(15) 



8(a,(3) 

is the Jacobian determinant for the transformation. 

Consider first a as defined by Vladimir. Suppose that 
f3 = 0; then by the Lorentz invariance properties of spectral 
radiance JO] we have 



Mfc) 
fc 3 



invariant 



(l + a)I(k) 
k 3 



and also 



m 
fc 3 



invariant 



(16) 
(17) 



if Vladimir and Estragon are to agree that the quantity 
I{k) admits interpretation as a radiance. The ratio 



(l + a)J(fc)_fi±^ 



Io(k) 



Io(k) 

is likewise invariant, so that 

1 + a = invariant, 



(18) 



(19) 



as it must be, as the ratio of two radiances evaluated in the 
same Lorentz frame. Thus 



a'(k') = a(k) = a(k') 
Next consider f3. We have 
(l + a)I + (3 



k 3 



= invariant 



from which 



i = invariant 
k 6 



(20) 



(21) 



(22) 



as must be for any radiance, in particular a noise radiance. 
By © we find 



P'(k') =/3'(r?fc) =rj 3 l3(k) 

The Jacobian is therefore 

~d{a',(3')- 



det 



L d(a,0) 

so Vladimir and Estragon must agree that 

g(a, (3)dad(3 = rfh{a' , f3')da'df3' 
and, by the principle of indifference PP,[S], that 
g(a 1 P))= V 3 g(a,r, 3 p)) 

with solution 

const 

g(a, (3)dad(3 — — - — d[3da 
P 



(23) 



(24) 



(25) 



(26) 



(27) 



Application of Bayes' theorem Q), as in 0, immediately 
gives the result that the joint posterior probability of T, e, 
a, and j3 is proportional to the product of @, J3J) , and 



P(T, e, a, (3, a \ I, K) oc exp 



(I -If) 2 



2<j 2 



, dT , d(3 dl 

d€ ^T da JV 
(28) 



V. Extended LST estimator 

Estimators for T and e may be constructed from H28J1 
as expectation values in exactly the same manner as in 1, . 
The treatment of spectral quantities integrated over a pass- 
band follows the equivalent discussion in £Q. In practice, 
will be unaltered for band-integrated radiances. 
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The calibration error parameters a and f3 are treated as 
nuisance parameters: One does not care what their actual 
values are, so long as they lie between specified limits. One 
is therfore at liberty to integrate (|28|l over those limits and 
obtain estimators for 



and 



(T) = 



J^-TP(T\I u a)^ 



and 



in terms of 



£Z eP((T),e| I it g)de 
P((T),e\I h a)de 



(29) 



(30) 



P(T,e\I,a) 



P(T\I,a) = 

err iicr p<7 i T > e > a ' & 



and 



(31) 



(32) 



In l|2"9")) and jSD) , (T) and (e) have no dependence on exactly 
what the calibration error parameters a and (3 were, for a 
given reported sensor aperture radiance. 

It does not appear feasible to integrate moments of (|2~5jl 
in closed form. However, by integrating over e first, it is 
possible to take advantage of the closed-form result for the 
LST posterior probability derived in P^: 



max i ^ nun 




) (33) 



(34) 



for each band i. As in pQ, 



k-2 



B k (T) 



i^(0)j ex P (-^)dk 



b = bib 2 



(35) 
(36) 



-Ft(0)ex P (-^)+ll(T,ii) )«//, 

7T jl 



(39) 



In (|33|l . as in l|28ll . the sensor radiance that appear in the 
quantites a, b, and c is related to the physical radiance by 

©■ 

The remaining integration over the nuisance variables a 
and (3 is now two-dimensional, and any integration over T 
to form the expectation value (T) makes for a third quadra- 
ture, for the full calculation. This is potentially awkward 
for routine evaluation, but the computational burden can 
be alleviated in special cases: 

1. If one knows a-priori that one is operating in a regime 
dominated by either gain or offset calibration error, the 
less important source of error may be ignored as a first 
approximation. 

2. In a vicarious calibration, the surface temperature 
may be accurately known. 

3. Once a value for (T) is obtained for one pixel in a 
dataset, expectation values (a) and (/3) can be calculated 
and used in estimation of (T) for other pixels. Should cal- 
ibration error be slowly varying, estimates of (a) and {(3) 
obtained from one dataset could be used for subsequent 
ones, or as initial guesses for updated estimates of (a) and 
(13). 
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with 



bi = 2 



fei 



B k (T) - -Fj:(0) 



exp{ 



)dk 



(37) 



fj, 



V k {r, l j)\dk-l, 



(38) 



